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Abstract. We present a complete polarization characterization of any quantum state of two 
orthogonal polarization modes, and give a systematic measurement procedure to collect the 
necessary data. Full characterization requires measurements of the photon number in both 
modes and linear optics. In the situation where only the photon-number difference can be 
determined, a limited but useful characterization is obtained. The characteristic Stokes moment 
profiles are given for several common quantum states. 
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1. Introduction 

Far from its source, any freely propagating electromagnetic field can be considered to a good 
approximation as a plane wave, with its electric field lying in a plane perpendicular to the 
direction of propagation. This simple observation is the root of the notion of polarization. 
At first glance, it may seem rather obvious how to translate such a concept into the realm of 
quantum optics. However, hurdles such as hidden polarization JT|, the fact that the Poincare 
sphere is too small to accommodate states with excitation larger than one photon |2), and 
the difficulties in defining polarization properties of two-photon entangled fields Q, to cite 
only a few examples, show that the classical theory, mainly based on first-order polarization 
moments, is insufficient for quantized fields. 

Here, we outline a systematic method for polarization characterization of quantum fields. 
The method is based on a simple premise; namely, that if we can predict the mth-order moment 
of the Stokes operator in any direction on the Poincare sphere, we know all there is to be 
known about the state polarization of this order Q, including any correlations between the 
Stokes operators. A tensor representation of the polarization information is based on such 
correlations. However, expressing the Stokes moments as functions of the measurement 
directions gives a more compact representation and provides a natural visualization. The 
Stokes profile representation also gives a relevant characterization for passive interferometry. 
Our analysis below makes use of both representations. 
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As a state's polarization properties do not require the full density matrix to be determined, 
it allows polarization tomography to be more easily performed than full quantum tomography 
|5]-[7]. Considering polarization tomography with ideal detection, we show that the number of 
measurement directions can be made equal to the number of independent parameters. 

The remaining material of the article is organized as follows. After recalling the 
fundamentals on the quantum description of polarization in section[2] we present our scheme 
for characterization of quantum polarization properties in section[3] In section@] we consider 
how the necessary data can be obtained experimentally. We thus arrive at an efficient method, 
which is feasible for polarization tomography of few-photon states. In section|5] we apply our 
characterization to several classes of states. Finally, our conclusions are presented in section|6] 

2. Setting the scene 

In the following, we consider monochromatic plane waves. Such fields can be decomposed 
into two orthogonal transverse modes, such as the horizontally and vertically polarized modes. 
For highly focused beams or waves in a waveguide, the plane-wave description is often 
inadequate, as the field is not longer transverse. It is our belief that the concepts discussed in 
this paper can also be extended to such non-plane waves, and several proposals have already 
appeared in the literature IBl fTOl . However, we shall not discuss such generalizations here. 

The classical theory for the polarization of plane waves was established by Stokes already 
more than 150 years ago IfTTl . We shall build on his theory as the basis of our treatment will 
be the Stokes operators, whose expectation values are the Stokes parameters 1121 . Following 
the conventions used in the quantum theory of angular momentum lfl3l and in quantum 
optics II14II15L we define the Stokes operators as 

(2.1) 

(, -jt at - \ o at - at n 

02 = lyanO-y — a^ay ) , 03 = a^au — a^ay , 

where a# and ay are the annihilation operators of the modes associated with horizontally and 
vertically oscillating fields, respectively. With this choice, the usual ordering of the Stokes 
parameters J* = (So), £ = (§3}, & = (§1), and Y = (§2) differs from that of the indices of 
the operators. However, as far as the theory below is concerned, we could just as well have 
associated any other pair of orthogonal polarization modes to these operators. That would 
only influence the interpretation of the theory and not the theory itself. 

As the annihilation and creation operators obey the bosonic commutation relations 
[a a ,a^] = 8 a p, for a,j3 € {H,V}, the Stokes operators satisfy the commutation relations 
of an su(2) algebra 

[Sj,S k ] = aejuS t , (2.2) 

where the latin indices run from 1 to 3 and Em is the fully antisymmetric Levi-Civita tensor. 
The noncommutability of these operators precludes the simultaneous exact measurement of 
the corresponding physical quantities. The variances (ASj) 2 = (5?) - (Sj) 2 are found to obey 
the uncertainty relation 

2(§o) < (A5!) 2 + (A5 2 ) 2 + (A53) 2 < (S (S + 2)). (2.3) 

Moreover, while the Stokes operators are all Hermitian, the noncommutability makes 
"mixed," non-symmetric products (such as S1S2) non-Hermitian, also precluding their direct 
measurement. 
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The standard definition of the degree of polarization for a quantum state p is 

where S = (Si, 02 j S3) and (S) is the Stokes vector. Note that only first-order moments of 
the Stokes operators are used in this definition. In a more elaborated characterization, the 
degree of polarization can be subdivided into excitation manifolds according to the total 
photon number N. This makes physical sense because since the corresponding observable 
So commutes with all the other Stokes operators 

M/]=0, (2.5) 

a complete set of simultaneous eigenstates of Sq and any of Si, §2, and £3 does exist. In 
fact, the statistics of the latter three operators is usually determined by a set of wave plates, 
a polarizing beam splitter, and two photodetectors, giving (in the ideal case) information not 
only about Si, §2, or S3, but simultaneously of Sq. 

Let us here take a quick look at excitation manifolds N = I and N = 2. One can readily 
convince oneself that any pure single-photon state j^i) satisfies Pj(|^Pi)) = 1, i.e., any such 
state is fully polarized according to the definition (12.41 . In fact, for an arbitrary single-photon 
state pi, the degree of polarization is related to the purity Tr(p 2 ) according to 



P S (P!) = v /2Tr( Pl 2 )-l. (2.6) 

However, this relation does not hold for other excitation manifolds. For example, any pure 
state in excitation manifold N = 2 of the form ifTBI 



|»F(a,e)) =ae- w \2,0) + iy/T^2cP-\l,\) +ae ie \0 1 2) , (2.7) 

where a and 9 are real numbers and < a < l/v2, satisfies P^d^a, 9))) = 0, which 
indicates that it is unpolarized. However, as we shall see below, these states have polarization 
structure (they are not isotropic in the polarization sense) and cannot be regarded as 
unpolarized. 



3. Higher-order polarization properties 

In order to characterize the polarization properties of a state, we shall employ measurements 
of higher-order moments of the Stokes operators. This is very close to the Glauber 
correlation functions in quantum coherence theory IfTTl , and has common grounds with 
Klyshko generalized coherence matrices JT|. In a recent paper {4), we have used the central 
moments for higher-order polarization characterization. Whereas the central moments may 
be preferred by some readers, the raw moments used in the present work seem to allow for an 
easier and more systematic approach. 

As we have already discussed, one can perform a measurement of the total photon 
number without disturbing the measurement of any other Stokes operator. In classical optics, 
this is tantamount to the fact that the state of polarization is independent of the intensity. This 
suggests that the polarization properties are given by S. However, an ideal measurement of 
polarization provides some information about the total energy and vice versa. For example, an 
even (odd) measured eigenvalue of any of the observables Si, §2, and S3 implies an even (odd) 
total number of photons. Also, determining the probability distribution for the total number 
of photons pn simultaneously sets bounds on the polarization properties in accordance with 
the inequalities ( 12.31 ). 
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Taking these observations into account, we distinguish polarization properties for 
different numbers of photons, and let full polarization characterization refer to complete 
knowledge of the expectation values of all possible combinations of the Stokes operators. 
The rth-order polarization information of a state p is then given by p^ and the expectation 
values of the form 

T hh-jr = $hh ■ -A)n = KifaSjA ■ - A). (3-D 
where jk € {1,2,3} and p^ denotes the normalized two-mode, Af -photon state obtained by 
projecting p onto the Mh excitation manifold 

P N =^. (3.2) 

Pn 

Using the Fock basis, the projector can thus be expressed as 1 # = L^=o \n,N — n){n,N — n\ 
and pn = Ty(Inp). For any given order r and excitation manifold N, the elements d3.ll ) form 
a Cartesian tensor T^'"^ (p) of rank r. Due to the Hermiticity of the Stokes operators, theses 
tensors satisfy 

T {r,N) _ r T {r,N) ,* n ^ 

h-jr 1 Jr-jl* ■ U ' ' 

We leave p^ and undefined for any N such that pn = 0, and employ the convention that 
they then do not contribute to sums. 

When (n,V,jk) is a cyclic permutation of (1,2,3), the commutation relation ( 12.21 ) 
implies that polarization tensor elements of neighboring ranks are related according to 

T (r,N) _ T (r,N) (r-\,N) 

h-W-jr-l h-VH-jr-l jl-jk-jr-1 ' ' ' 

Hence, Iv -1 '^ can be determined from and, consequently, T^ ,Ar ) determines all TWO 
such that r < R. Complete polarization information of order R is thus equivalent to complete 
polarization information of all orders r <R. 

Using the relations d2.U and ( 12.2b . it is also straightforward to show that the polarization 
information carried by T( R - N ) is equivalent to that contained in the set of generalized coherence 
matrices of orders 2r {r < R), whose elements are of the form [Q~) ((a#) J ' (a v ) r ~ ' ' a k H a r y k ) n . 
Having complete polarization information of all orders about a state p is therefore equivalent 
to knowing its block-diagonal projection 



Ppol = £ PNpN, (3-5) 

N=0 

where pV is given by ( 13.2b . This is the so called the polarization sector (or polarization 
density matrix). The number of parameters characterizing a block-diagonal state limited to 
the excitation manifolds Ni , N2 , ■ ■ ■ , N v is 

-1 + £(A^ + 1) 2 . (3.6) 

k=] 

In particular, when a state is limited to the manifolds 0, 1, . . . ,N, the number of parameters 
simplifies to 

-l + I(A- + D 2 = ^ 2+ f +13) - (") 
Co 6 

For such a state, complete polarization information of order TV is sufficient to determine its 
block-diagonal projection ( 13.51 ). The general density matrix for a state with no more than N 
photons is determined by N(N + 3)(N 2 + 3A^ + 4)/4 independent real numbers. Hence, the 
polarization share of this information quickly decreases with N. 
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4. Polarization tomography 

We now turn to the question of how to characterize polarization properties experimentally. 
Since our interest is limited to the information contained in the block-diagonal 
projection ( 13.51 ). it is clear that we are not required to do full quantum state tomography j5]-[7). 
As complete polarization information corresponds to doing quantum tomography of all N- 
photon Hilbert spaces excited by the considered state, one can make use of the methods 
developed for finite-dimensional systems |fT9H2Tl . However, some recently proposed higher- 
order intensity measurements El seem to be closest related to the ones we present below. 

We will assume ideal measurements and that the total photon number and its probability 
distribution p^ can be determined. This is obviously a severe restriction apart from the lowest 
excitation manifolds. However, the situation where no information about the total photon 
number can be obtained is described by simply summing over the different manifolds as 
discussed in section R~4l 

Below, we also treat the experimental determination of different moments of an 
observable as different measurements. In principle, each moment requires an infinite number 
of measurement runs in order to be determined exactly. This would also give us the full 
probability distribution of the eigenvalues and thus all the moments. However, for the vast 
majority of realistic probability distributions, a lower moment requires fewer runs to be 
accurately determined. 



4.1. Moment measurements 

The fact that the classical Stokes parameters are easily determined experimentally makes 
them highly practical. Also in quantum optics, the measurement setups corresponding to the 
fundamental Stokes operators are simple. These setups are composed only by phase shifters, 
beam splitters and photon-number measurements. The effects of linear optical devices are 
described by SU(2) transformations fl4l . which can be expressed as 



#(3>,®,E)=e 



-;$S3/2 e -,-©52/2 e -;ss 3 /2 



(4.1) 



where <t>, ©, and E are the Euler angles. Any such transformation can be easily realized using 
linear optics |23l and they are lossless, so they leave Sq unaffected. 

Let us now introduce the Stokes operator in an arbitrary direction characterized by the 
unit vector n £ I 3 as 



S a = n S 



k=l 



The effect of an arbitrary SU(2) transformation on S n , can then be expressed as 



U{<S>,®,Z)S n U\<S>,®,Z)=S-. 



R3(*)-R 2 (e)-R 3 (H)-ii) 



(4.2) 



(4.3) 



where R<.(0) denotes the matrix describing a rotation of <p around the e^-axis, e.g. 

"10 

Ri(</>)= cos</> -sin</> . (4.4) 
sin <p cos 

Hence, any SU(2) transformation corresponds to a proper rotation in M 3 fl4| . We note that 
£3, which gives the photon-number difference, is transformed according to 

&(<Z>,®,S)S 3 f (<Z>,®,Z)=S n , n= (sin0cos<I>,sin0sin4>,cos0). (4.5) 
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That is, U(<P,&,0) is the standard displacement on the sphere and the transformation 
parameters and <t> equal the spherical coordinates of the vector n characterizing the 
transformed Stokes operator. We see that any S n is related to £3 by an SU(2) transformation 
corresponding to a polarization rotation of 0/2 followed by a differential phase shift of <E>/2. 
Expectation values of the form (S' n ) can thus be straightforwardly determined experimentally. 
Ideally, we can simultaneously measure the total photon number N, so that also expectation 
values of the form (S'„)n can be determined. 

The tensor gives any expectation value of the form 

(Sn l S B2 ...S nr ) N = t ...Inf.-n^S. (4.6) 

71=1 j r =l 

When all vectors are the same, ( 14.61 ) simplifies considerably. For a given state, the relation 
between (S' n ) n and the direction n will be referred to as the -photon Stokes moment profile 
of order r. These profiles can be expressed as 

(s' n )N = i L'^K k Mf ] , (4.7) 

k=oe=o 

where the moment component M^f^ is the sum of all tensor elements of the form T^' N j r 
that have k ones and £ twos as subscripts. Due to ( 13. IK every moment component is thus 
the expectation value of the Hermitian operator formed by the sum of the Stokes-operator 
products corresponding to its tensor elements. The number of such elements is given by the 
trinomial coefficient 

(k,£,r-k-£)\ = — - -. (4.8) 

v ' ' ' kW.(r-k-£)\ 

For example, we have = T§f> + + T gf> + T§?> + + . We note 

that the sum of all tensor elements of order r in excitation manifold N can be written as 

t-t T h"l = t L< ] - 3 r/2 <^> . (4-9) 

/l=l j r =l k=0t=0 

where n diag = (1,1, l)/v / 3. 

Since the polarization tensor satisfies the Hermiticity condition (13.3b . the moment 
components are real, and consequently the real part of is sufficient to determine n in 
any direction n. Naturally, knowing the Stokes moment profile ( 14.7b is equivalent to knowing 
the 

(r+l)(r + 2) 

m r = (4.10) 

moment components M^f' . 

Moreover, using the commutation relation (12.2b . it is possible to determine the 
differences between the tensor elements belonging to the same moment component MXf' 
from Tv~ l ' N ). Since every element of belongs to such a moment component, it thus 
follows that T^'"- 1 ^ together with all M { k r f ] determine T*''^. 

Let us now introduce a standard ordering of the Stokes operators according to 

0^J=S\s{S r f k - 1 . (4.11) 

Making repeated use of the commutation relation ( 12.2b . the moment components can then be 
expressed as 

= (k,£,r-k-£)\ {&Q) N + {<fj$) N , (4. 12) 
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where is a sum over Stokes-operator products of orders smaller than r. The Casimir 
operator 

3 

S 2 = £^=S (S + 2) (4.13) 

k=\ 

implies that, for r > 2, we have 

<^+2> + {&$ +2 )s + {^})s = N(N + 2){^- 2) ) N + {S\ [^,si]S^- f - 2 ) N , (4.14) 

where the last term can again be written as a sum over Stokes-operator products of orders 
smaller than r. Equations (I4.12l i and ( 14.14b show that there is a relation between moment 
components of orders r and r — 2, and that the number of independent moment components 
of order r is m r — m r ~ 2 = 2r + 1. For a general TV-photon state, the number of independent 
moment components to determine is thus L^Li(2r+ 1) = N(N + 2). For a general block- 
diagonal state, we also have to determine the probability distribution for the total number 
of photons. Assuming that the excited manifolds are known to be limited to N\ ,N 2 , • • • ,N V , 
we find the number of independent parameters to be v — 1 + YX=i^k(Nk + 2), which is in 
agreement with (13 .6b - 



4.1.1. General single-photon state In the basis (| 1,0), |0, 1)), the density matrix of a general 
single-photon state can be written as 

where R 2 +I 2 < jIq(1 — Jib). Using a superscript to identify state-specific average values, the 
first-order Stokes moment profile is given by 

(£„)?' =27?rti-2/n 2 + (27io-l)«3- (4.16) 
Hence, the three moment components are seen to be independent. 



4.7.2. General two-photon state Using the basis (|2,0), |1,1), |0,2)), the density matrix of a 
general two-photon state can be written as 

/ Jii Ri + ih Ri + ih \ 
p 2 = Ri-ih Ji 2 R3 + U3 ■ (4.17) 
\ R 2 - ih R3 - ih 1 - JTi - JT 2 / 

The first- and second-order Stokes moment profiles can then be expressed as 

(S a f 2 2 = 2V2 +R 3 )m - (h +h)n 2 ] +2(2^i +7C 2 - l)n 3 , (4.18) 



{Sl) p 2 2 = 2 (1 + ji 2 + 2R 2 )n\ +2 (1 + 7t 2 - 2R 2 )nj + 4(1- 7t 2 )4 - 8I 2 nm 2 
+ 4V2n 3 [(i?i-i?3)ni - [h-h)n 2 ], (4.19) 

which makes it easy to identify the moment components. As implied by ( 14.131 ), the moment 
components of the three first terms of ( 14.19b are determined by two parameters. Hence, there 
are only five independent second-order moment components. 
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4.2. Choosing measurement directions 

We have seen that the information content of the moment components allows us to do 
polarization tomography by only measuring moments. Performing the moment measurements 
in increasing order, the 2r+l independent moment components for each order r and manifold 
N can be determined by choosing equally many directions n such that ( 14.7b gives linearly 
independent equations for the unknown moment components. 



4.2.1. First order Quite naturally, both the first-order moment components and the first- 
order polarization tensors are given by the manifold-specific Stokes parameters 

M^l=T^ = {§j) N . (4.20) 

The three sets of information are hence identical and are obtained by determining the 
expectation value (S n )N for the directions n = (1,0,0), (0, 1,0), and (0,0, 1). As the operators 
S n and S- n only differ by the signs of their eigenvalues, the corresponding measurements will 
give the same information. Hence, equivalent measurements correspond to a line through the 
origin. Choosing three orthogonal directions as above thus results in a uniform distribution of 
the measurements on the Poincare sphere. 



4.2.2. Second order We have seen that there are five independent moment components of 
second order. Thinking of the measurements as lines, we choose the directions as 



ni,2 : 



(0,±2,1 + V5) 



1*3,4 



(±2,1+ Ao) 



n 5 



(1+ Ao,2) 



(4.21) 



^10 + 2^ ^10 + 2^ \^W+2V5 

which maximizes the minimum angle between the lines l24ll25ll and thus in some sense 
spreads out the measurements over the Poincare sphere as much as possible. The six second- 
order moment components are then given by 



M, 



(2,A0 



0.1 



V5 (§ 2 

2 { n ' 



J B 2 I 



M 



M 



(2,A0 
1,1 

(2,A0 
1,0 



V5 
2 
V5 



(si 



-sl + sl 



M, 



(2.A0 
0,0 



(15+7A)(5i+5i 



VSN(N + 2) 
(lQ + 4V5)(S 2 n3 +S 2 a ^ 



{6 + 2V5)N(N + 2) 



M, 



4(7 + 3 y/5) 

_ (10 + 4V5) (§ 2 ni + Sj 2 ) N + (25 + 11 A (§ 2 nj + § 2 n4 ) N - (14 + 6 V5)N(N + 2) 



0.2 



M {2 ' N) - 

M 2,Q — 



4(7 + 3 V5) 

(36 + 16 V5)N(N + 2)-(25 + nV5)(S 2 ni +S 2 n2 ) N -(l5+lV5) (§l 3 + §i 4 ) 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



,(4.26) 



-(4.27) 



4(7 + 3 y/5) 

Note that their determination does not require any first-order measurement. The second-order 
polarization tensors can be expressed in the first- and second-order moment components as 



T (2,A0 



Ml 



(2,W) 



M 



(2jV) 



A/ 



(2.,V) 



2.0 



■ + i(§3) N 



M 



(2,W) 



2 

(2.JV) 



~ i{S3}N 
+ i(S~2)N 



M 



M, 

(2,N) 



(2-N) 



2 

(2jV) 



-i(S: 



\ 



2)N 



0.2 



i(S~i)N 



M, 



(2,N) 



(4.28) 



0.0 
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When writing tensors, we let larger entities and rows correspond to tensor indices placed to 
the left of those corresponding to smaller entities and columns. 

As an aside, we decompose the Stokes-operator covariance matrix into different 
excitation manifolds T = Y*n=oPnTn and note that the matrix elements are given by 

_ (SjS k ) N + (S k Sj) N ~ ~ -U P (T<< 2 - N h T (l,N) T (l,N) 9Q s 

Ljk.N= ^ {!>j)N{!> k ) N -Re[l jk l k . (4.29) 

For any state that satisfies P$ = 0, we thus have = Re(T^ 2 ' N ^). 



4.2.3. Third order We know that there are seven independent third-order moment 
components. Maximizing the minimum angle between seven lines 11251 . we find that the 
measurements should correspond to S\, §2, S3, and the directions 

(±1,1,1) _ (±1,-1,1) 



114,5 = 



V3 ' ^~ V3 ■ 

However, this choice gives only four independent measurements, since we have 



S 



3y/3 



(-SL + 



3N(N + 2)-4 



(4.30) 



(4.31) 



and similar relations for S\ and S\. By choosing directions close to S\, §2, and S3, it is possible 
to determine all third-order moment components. However, this choice would make it hard to 
obtain the necessary data, since the corresponding expectation values differ only slightly from 
the known (S\)n, (S^)n, and (S^n- Consequently, although highly symmetric polyhedrons 
have been successfully applied in protocols for tomography of multi-qubit states [26], the 
related method considered here fails. How to optimally choose the measurement directions 
for higher-order polarization tomography thus appears to be a complicated problem. This 
notwithstanding, the third-order polarization tensors can be expressed as 



T (3,A0 = 



M 



(3,W) 



(2JV) 
3,1 



{3,N) (2,N)_ 
2.0 1,2 



M.N) 



J3.N) 



,(3,«) 



■+iTj 

(3,W) 



(2,«)_ (2.JV) (2,W) 



u 2,2 +,J 3,3 




3" 

2.JV)_ L(T (2,«) 



J3.N) 



I.I 

+i4T. 



(2.N) 



2,2 
+i2K 



" !T 3, 
(2,W) 

2J 



(2./V) 



— er, 



(2.W) +(T (2.JV)^. T (2,,V) 



I.I 



2.2 



(3JV) (2JV) (2,/V) 
M 1.0 "'3,2 +,/ '2,3 



J3.N)_ 



7F- 



3.3 



1,2 

J2.N) 



,(3,iV) 



+;2r. 



cm 
,1 

3 — 

(3JV) 



.(4.32) 



4.3. Recurrence relation for Stokes moment profiles 

Above, we have seen that the N lowest-order Stokes moment profiles {(S^)n}^ = i contain all 
polarization information of an -photon state. In particular, we show in the appendix that the 
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higher-order profiles are determined by the recurrence relation 



(5™ 



JV+l+jU\ 



N/2 



£ 4 N / 2+i -Jf(N + 2,2j) (Sl j - 1+ ») N , N even, 



= < 



7=1 



(4.33) 



-^4^-Jf(N + 2,2j+l)(S 2 n 



/Vodd, 



where jj. is a non-negative integer and f(n,k) are the central factorial numbers of the first kind 
given by 1271 

0, n < k, 

<5„,o, k = 0, 

n-k i ,v; /->,„ ; / ;\ / ; \ n-k+j (4.34) 



2n-A "p* (-1V (In-lk 



/(«,*) = 

f~bj l ( n + j) \n-k-j 
For the lowest excitation manifolds, we thus get 

(%>o = 0, 



K-ir 

»!=() 



, !<£<«. 



2 r - 2 (52) 2 , reven, 
2 r - 1 (5„) 2 , rodd, 



(S'n) 



1, reven, 
(5„)i, rodd, 



9-3 r +(3 r -l)(^) 3 
8 



(4.35) 
r even, 



(9-3-M(A 1 )3+(3-'-l)(^) 3 rodd _ 



The above relations for N = I and = 2 can also be established from the general property 



5jn, reven, 
— rodd. 



(4.36) 



4.4. Non-resolved photon numbers 

As pointed out above, apart from the few lowest excitation manifolds, it is difficult to 
distinguish different excitation manifolds experimentally. In case there is no information 
about the total photon number available, the measured expectation values are weighted 
averages over the manifolds of the form (A) = Y.n=qPn(A)n- Due to linearity, it is clear that 
the corresponding polarization tensors, Stokes moment profiles, and moment components, 
which are given by 

TW = £ p N T^\ (S' n ) = £ p N (S r a ) N , Mg = £ p N M^\ (4.37) 

N=Q N=0 N=0 

respectively, enjoy most of the properties of their manifold-specific counterparts. However, 
due to the factor N(N + 2) appearing in ( 14.14b . the relations between moment components of 
different orders depend on the excitation manifold. In general, this makes all photon-number- 

(r) 

averaged moment components M k I independent and /?th-order polarization tomography then 
requires \ m r = R (R 2 +6/?+ll)/6 parameters to be determined. However, the knowledge 
of the average photon number (So) and its variance (Sq) — (So) 2 is sufficient to remove the 
redundancy for the second order, since we then know the right-hand side of the relation 
M^l +M^l + M ( f ( j = (Sq(Sq +2)) obtained from j4~T3l . With this partial knowledge about 
the photon distribution, we can thus determine the second-order moment components using 
the five measurement settings given in section |4.2.2| 

Now, assume that we know that a state is limited to the first three manifolds, i.e., that 
the number of photons cannot exceed two. In this case, the determination of (Sq), (Sq), 
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and the three lowest-order Stokes moment profiles is sufficient for complete photon-resolved 
polarization characterization. Explicitly, we have p\ = 2 (So) — (Sq) and pi = ((So) — (§o))/2, 
which together with (14.35b give 

,*v 4&) - (§1) * (§1) - (5n) «. 2((S 2 n ) + (Sp-2(S )) 

^ a ' i = 77F~\ — T7&\> \ ! >nh = = , VV2 = t^t — 77T7 ■ ( 4J8 ) 

6(S )-3(Sg) 3 (So)-( s o> 



5. A menagerie of states and their polarization properties 

We next apply the characterization developed above to some classes of states. In most cases, 
we give only the Stokes moment profiles for the states, as these provide the most compact 
presentation of the polarization properties. It should be straightforward to obtain the moment 
components and polarization tensors if needed. 

Using ( 14.5) . we can write the Stokes moment profiles of an arbitrary state p as 

(§' n ) N = Tr[J/(4>, 0, S) § r 3 0t 0, S )p], (5.1) 

where the SU(2) transformation t/(<t>,0,E) is given by (14.1b . Now, consider the state p' 
obtained by applying an SU(2) transformation to p according to 

p 1 = U{<p,V,t;)pu\<p : V,t;). (5.2) 

As the trace of a product is invariant under cyclic permutations, (14.3b ensures that the Stokes 
moment profiles of the state p' are related to those of p by rotations. Indeed, we find that 

(S' n )^ is obtained from (SQ^ by rotating the latter £, around the e3-axis, followed by a rotation 
of # around the e2-axis and another of <p around the e3-axis. That is, we have 

{Sn) P N = (^R 3 (-5)-R 2 (-d)-R 3 (-(p)-ii)Af ( 5 - 3 ) 

Naturally, the sequence of rotations appearing in ( 15.31 ) is the inverse of the one described 
above. Because of these simple rotations, the determination of the polarization properties 
of a state p, implicitly gives the polarization properties of all states related to p by an 
SU(2) transformation, although these states may appear very different. We note that 
common, passive, two-mode interferometers are described by SU(2) transformations too. The 
considerations below are therefore relevant to interferometry. 



5.1. SU(2) coherent states 

The SU(2) coherent states are the eigenstates of the operators S n . They are also the only 
states that minimize the variance sum, i.e., that saturate the left inequality in the uncertainty 
relation ( 12.3b . Using the spherical coordinates (14.5b , we have the eigenequation5 n |A^; ©,<!>) = 
Af|Af;©,4>). The jV-photon, SU(2) coherent states are of the form 



N 



=0 




|AT;0,4>) = £ e~^ x j ( ^ ) sin*- ( - ) cos" ( - ) \n,N-n) 







= e- im / 2 U{<P,e,0)\N,0). (5.4) 

Since an overall phase factor does not have any physical significance, they are thus all related 
to the state \N, 0) by an SU(2) transformation (14.1b . As discussed above, such transformations 
correspond to simple rotations of the Poincare sphere, so we limit our explicit treatment to 
the states \N, 0) . Making use of well-known results for the beam splitter, we easily find the 
Stokes moment profiles to be 

(S^ = I(*-2*)'(T) sin- ( f) cos^-> ( f) . (5.5) 



k=0 
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The lowest-order polarization tensors are 





T^flJV.O)) 










iN 2 


N(N-2) 


iN(N-2) 











-iN 2 








N 2 


-iN(N-2) 


N(N-2) 





N 2 


iN 2 





-,7V 2 


N 1 






\ w 3 / 

For Af > 0, the SU(2) coherent states thus satisfy P s = 1. 



N 


/W 





-iN 













yv 2 



(5.6) 



(5.7) 



5.2. Two-mode coherent states 

Since any pair of two-mode coherent states with the same average total energy are related by 
an SU(2) transformation, it suffices to study states of the form 



l«,0> 



N=Q 



a" 



/AH 



\N,0). 



(5.8) 



The block-diagonal projection is clearly independent of the phase of a, and is given by a 
Poissonian mixture of SU(2) coherent states that all belong to different excitation manifolds. 
Hence, the manifold-specific expectation values coincide with those of the SU(2) coherent 
states, and the manifold-averaged Stokes momentum profiles ( 14.371 ) become 

«) |a ' o> =£^^>r > , (5.9) 

N=0 



N\ 



where N = \oc\ 2 . Using the corresponding tensor relation (14.371 and results for the states 
|Af,0), we easily obtain 



(Sn, 



T«(|a,0» 



l«,o) _ 



ATrt3 , 



02^.0) =N ( l+Nn 2^ 





I , T( 2 >(|a,0)H 
A' 



N 

-iN N 



(§3}\afi) = N n3 ^i + 3N+N 2 nl), (5.10) 
iN 

(5.11) 







T (3) (|«,0)) 





N(N+1) 



N(N+1) 
\ 









iN(N+l) 


N(N-l) 


iN(N-l) 











-iN(N+l) 








N(N+l) 


-iN(N-l) 


N(N-l) 





N(N+l) 


iN{N+l) 





-iN(N+l) 


N(N+l) 






(5.12) 



\ N(N 2 +3N+l) ) 

In accordance with classical optics, = 1 for any two-mode coherent state with a finite 
average photon number N. We note that when N ^> 1, the lowest Stokes moments satisfy 
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(5J,)' a '°^ | j w (Nn^Y . For «3 = 1, the Stokes moments ( 15.10b are directly given by 
the Poissonian photon distribution, and the approximation corresponds to the classical 
deterministic limit. The «3 dependence of the approximation describes the transmission 
through a classical beam splitter or linear polarizer. In particular, Malus' law is obtained 
for r = 1 . 



5.3. \m,m) states 

When considering the two-mode Fock states \m,m), which allow for Heisenberg- 

limited interferometry l28l . we implicitly treat all states obtained from these by SU(2) 
transformations. The latter states can be expressed as 

2m 



fit a e\i v sin'"z? £ 2V(2m-t)!t! 
U(<p,-&,£) \m,m) = — — — > — — — ; — 



m\2 m ^ tan^tf 
V S}(m\( m-j \ / tan 2 #V' f /e t „ 

where [x\ denotes the largest integer that is smaller than or equal to x. We have also assumed 
that the binomial coefficients are defined through the gamma function, so that negative 
integers are allowed as arguments. The Stokes moment profiles of the states \m,m) are given 
by 

= { ^D(2;-l)!!]V(,2;)( m 2 y')sin^ 0I .even, (f . M) 

0, r odd, 

where F(n,k) denotes the central factorial numbers of the second kind iA.2i . As both 
arguments are even in our case, we have 11271 

, . J, (~iy +k k r 

y( ^- 2 E (A*w-*)i - (5 " 15) 

In particular, we get 

<©]P-5«*±2** ,5.16) 



N n 



> _ ,rt,r , ^ • 2« 16 + 3(JV-2)(jV + 4)sin 2 © 



(Sj)]7' 7/ = N(N + 2)sm z ® i 1 . (5.17) 

The photon-number symmetry makes all odd-order Stokes profiles vanish. This is in stark 
contrast to the states |iV, 0) and \a,0), and makes P$ = 0. However, the so-called hidden 
polarization of the \m, m) states appears in the even-order Stokes profiles. As sin 2 © = 1 — « 2 , 
these profiles have common features with rotated ones for \2m,0). This is in agreement with 
the fact that more elaborate measurements than those considered here are required to achieve 
Heisenberg resolution when employing the \m,m) states 



5.4. Two-mode squeezed vacuum 

Using the process of spontaneous parametric down-conversion, one can straightforwardly 
generate two-mode squeezed vacuum states. These have thermal photon-pair distributions 
and take the form 

00 / 2N'" 

l*TM S v)=I>*y {2+ff)m+i Km), (5.18) 
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N = 2 




Figure 1. The Stokes moment profiles (5„ ) N for the W-photon NOON states with 1 < N < 6. 
Dark and light surfaces indicate the directions for which (S^)n takes positive and negative 
values, respectively. 



where N denotes the average number of photons. From ( 14.371 ). ( 15.161 ) and ( 15.17b . we thus get 

( ^TM S y = iV(2iV+3)sin 2 @; ^ 

^)^sv =j v sin2Q 32 j V + 48 + 9(8 j y3 + 20^ + 10^-5)s 1 n^_ ^ 



5.5. NOON states 

Finally, let us consider the NOON states (\N,0) + \0,N))/y/2, where N > 0. These are in 
some sense optimal for interferometry l30l . Their Stokes moment profiles are found to be 

r o, 



(S r n 



r\NOON 
N 



4(W-l)/2 



k=0 



r odd, N even, 
r and N odd, 



k=0 



N 



„2k 







1) cos (NO) cos [ — | 







Jl{N-k) 











(5.21) 



We note that when r is even and N is odd, the effect of the second term within square brackets 
vanishes. For each of the first six NOON states, we have plotted the Stokes moment profile 
(5„)jv in figure [TJ In the horizontal plane «3 =0, we have 



W\NOON| 



/cW\JN' 



Nlcos(N®)+2 N / 2 Q N/2 (N/2), N even, 



I &=7t/2 



(5.22) 



N\cos(N&), 



N odd, 
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where the polynomial 1311 

Qj(n) = 2J- 2nl £^ k j(n-k) 2 J (5.23) 

satisfies Qo(n) = 1 and the recurrence relation Qj+\(n) = 2n 2 Qj(n) — n(2n — l)Qj(n — 1). 
Hence, the de Broglie wavelength of the NOON states, which scales as N~ 1 and is the reason 
for their superiority |f32l , can be seen in these measurements. 

Let us now return to the pure two-photon states that satisfy F$ = 0. These are given 
in (12. 7t and are found to be related to the two-photon NOON state by appropriate SU(2) 
transformations according to 

|<PM)) = ?/(| + 1 arccos(V2a),-|) (5 .24) 

Hence, any Stokes moment profile of the state ^(a,©)) is related to the corresponding 
NOON profile by simple rotations. Since R 3 ( n/2 + 6 ) • R2U) • R3 ( - n/2) = R 3 ( 6 ) • Ri {-%) , 
(^b/z" i s obtained from (S^)^ OON in figure Q]by applying a rotation of — arccos(v / 2a) 
around ei followed by a rotation of 6 around 63. Since (S n ) 2 = is independent of 
n, the states ^(a, 9)) lack first-order polarization structure. However, they do all have a 
second-order polarization structure. 



6. Conclusions 



Using expectation values of Stokes-operator products, we have developed a systematic 
scheme for characterizing higher-order polarization properties of two-mode quantized fields. 
Polarization tensors and Stokes moment profiles were introduced as two representations of 
the polarization information. The latter show how passive interferometry affects the moments 
of photon difference. This viewpoint was taken as polarization properties of different states 
were compared. 

Other possible representations of the polarization information include central moments 
H, quasi-probability distributions l33l and excitation-specific generalized coherence 
matrices. Complete polarization characterization requires the excitation manifolds to be 
addressed separately. For situations where this cannot be achieved, our characterization 
coincide with the one provided by Klyshko's generalized coherence matrices 0]. 

Assuming ideal photon-number resolving detectors, we have shown that it is possible to 
efficiently collect the data through Stokes moment measurements in different directions. In an 
experiment, it may be more practical to use more measurement directions than the minimum 
required, but our method should serve as a guide. In particular, we expect the introduced 
moment components to be useful. Another advantage of the described method is that it treats 
the Stokes moments order by order. Hence, if only the first few polarization orders are of 
interest, it makes the measurements easier. 

Since the different excitation manifolds are treated separately, losses have drastic 
consequences in that higher excitation manifolds then contribute to the lower ones. 
Furthermore, whereas linear losses often model imperfections of single-photon detectors well, 
photon-number resolving detectors, which are required for full polarization characterization, 
are more complex and may call for nonlinear modeling. 

On the other hand, the separation of data into excitation manifolds and moment 
orders may be useful when developing methods for efficient determination of polarization 
characteristics. For example, one can take into account that all state projections in the 
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different excitation manifolds must be physical states. In this way, it should be possible 
to develop efficient maximum likelihood methods similar to those regularly employed in 
common quantum tomography. 
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Appendix 



This appendix gives a derivation of the recurrence relation ( 14.331 ), which involves central 
factorial numbers |27l . We let the sets of non-negative and positive integers be denoted as No 
and N+, respectively. For x <G R, the central factorial of degree n is defined by 

1, n = 0, 



n-2 

* n 

k=2-n 



«eN+. 



(A.l) 



The central factorial numbers of the first and second kind, f(n,k) and F(n,k) with n,k G No, 
respectively, are then defined through the expansions 



jcM = £/(n,*)jc* 

k=0 



x" = £F(n,*)*M. 

A-=0 



We note that 

f(n,0)=F(n,0) = S nfi , 

/(«,«) =F(n,n) = 1, «eN+, 

where 5„i denotes the Kronecker delta. For x <G R, we clearly have 

v-1 



r [2v] 



7=0 



V G 



x [2v+l] 



7=1 



X -l J "2 



V € Nq. 



(A.2) 



(A.3) 
(A.4) 



(A.5) 



(A.6) 



Consequently, both f(n,k) and F(n,k) vanish if one argument is even and the other is odd. 
For m e N , V G N+, and m < V, (IA 2b . dA.3b and (IA.51 > give 



,Pv] 



Factoring out an m and making use of ( 1A.4| |. we obtain 



^/(2v,2;)m 2/ = 0. 



2v-l 



v-l 



■£/(2v,2j)m 2M - 

7=1 



(A.7) 



(A.8) 
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Similarly, for m, v € N+ and m < V, it follows from (IA.6t that 

(m - 1 /2)P V+1 ] = £ /(2v + l,2j + l)(m- 1 /2) 2 ^' +1 = (A.9) 

and 

£/(2v + l,2j + l)(m-l/2) 2 > = 0. (A. 10) 

.7=0 

Hence, for any integer m satisfying 1 — V < m < V, we have 

v-1 

(m - 1 /2) 2v = - £ /(2v + 1,2; + l)(m - 1 /2) 2 >. (A. 1 1) 

;'=0 

Now, consider an observable A in <f -dimensional Hilbert space with eigenvalues {A 7 }^ =1 . 

In the eigenbasis, we then have A r = Diag(A[, A£ , . . . , AJ). If all eigenvalues are integers and 
satisfy < V, where v £ N+, (I A.81 > gives the recurrence relation 

v-1 

pv-i +M = £ /(2v,2j)A 2 ^ 1 +' 1 , (A.12) 

7=1 

which is valid for any fi <G No- If all eigenvalues are half-integers and satisfy |A,t| < V — 1/2, 
where v e N+, it follows from Eq (IA.1U that 

v-1 

A 2v+ V = -Y J f{2v + l,2j + l)A 2j+ ». (A. 13) 

j=o 

Here, fi can take any integer value, since the inverse of A is guaranteed to exist. We note that 
we can apply arbitrary unitary transformations to both sides of (|A.12t and (|A.13t , so they are 
valid in any basis. 

Now, consider the angular momentum component Jt, of a spin-.? particle. The dimension 
of the corresponding Hilbert space is d — 2s + 1 and the eigenvalues are — s, — s + 1, . . . , s. For 
s £ N+, the magnitudes of all eigenvalues are integers less than s + 1 and (IA. 121 > becomes 

Jf +1+ " = -if(2s + 2,2j)jT 1+ ^ (A.14) 
7=1 

If s is not an integer, all eigenvalues are half-integers with magnitudes less than or equal to s 
and (IA. 1 31 > becomes 

4 S+1+ " = -~ji 2 f(2s + 2,2j+l)% + ». (A.15) 
;'=o 

Taking the expectation value of the corresponding recurrence relations for the Stokes operator 
£3 = 2/3 in excitation manifold N = 2.?, gives the desired result (I4.33l l. 
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